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Abstract 

Plasma instabilities are parametrically the dominant nonequilibrium dynamics of a weakly cou- 
pled quark-gluon plasma. In recent years the time evolution of the corresponding collective colour 
fields has been studied in stationary anisotropic situations. Here I report on recent numerical 
results on the time evolution of the most unstable modes in a longitudinally expanding plasma 
as they grow from small rapidity fluctuations to amplitudes where non-Abelian self-interactions 
become important. 



The experimental findings at the Relativistic Heavy Ion Collider RHIC [T] suggest a 
much faster thermalization than can be accounted for by the (weak-coupling) bottom-up 
thermalization scenario of Ref. [2]. However, as first pointed out by Arnold et al. [3], this 
scenario did not include the inevitable presence of plasma instabilities in an anisotropic 
quark-gluon plasma |4|5j . In fact, these are even the most important collective effects 
to leading order. Much effort has by now been invested recently in studying the specific 
non-Abelian dynamics of plasma instabilities, in particular in the crucial regime where 
the associated colour fields are still too weak to affect the distribution of hard particles 
but already strong enough for essential non-Abelian self-interactions. 

Initial numerical studies that concentrated on the most unstable modes, namely long- 
wavelength coherent colour fields that are constant in transverse directions, suggested 
continued exponential growth also in this nonlinear regime [5J. Subsequent studies [718] 
with fully 3+1-dimensional colour fluctuations found however a late-time transition from 
exponential to (rather weak) linear growth. On the other hand, the more recent simula- 
tions of Ref. [9] found a continued exponential growth of initially small perturbations in 
the case of very strong momentum anisotropy, at least for small initial field configura- 
tions, which is similar to the behaviour of the 1+1-dimensional case. 

All the above-mentioned investigations were performed within the hard-loop effective 
theory of a collisionless stationary anisotropic plasma [10] . Clearly, in view of the crucial 
dependence of the fate of non-Abelian plasma instabilities on the amount of momentum 
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anisotropy, which rapidly increases during the free-streaming stage, a generalization to 
Bjorken expansion is desirable PI 

In the stationary anisotropic situation, the coloured fluctuations of a colour-neutral 
space-time-independent background of hard particles can be factored according to 

Sf a (x;p) = -gWZ(t,x;v)dfcf (p) (1) 

where the auxiliary fields W depend only on velocities tj m = p M /|p| = (l,v) and are 
governed by a Vlasov equation 

[v ■ D(A)]W^x;w) =F m {A)v^ (2) 

in which the hard momentum scale p° does not appear. The latter can be "integrated 
out" , determining just the mass parameter in the induced current j M [6f] — m 2 D J v W(x; v) , 
where W is one particular linear combination of the components W u . The collective colour 
fields associated with plasma instabilities can then be determined by solving the equation 
for W together with the non-Abelian Maxwell equations D p (A)F ptl — j^{x). 

Remarkably, this formalism can largely be taken over to the nonstationary case of a 
plasma undergoing Bjorken expansion |12j . 

In comoving (rapidity) variables x a = (t = \/t 2 — z 1 ,x l ,rj = atanh(z/i)) and p a = 
|p_i_|(cosh(y — 77), cos sin </), T _1 sinh(y — 77) with y = atanh(p°/p 2 ), a free-streaming 
background distribution is given by 



/o(p, x) = /iso {\Jp\_ +P 2 n lrl^j = / lso {\Jp\ + (p' z r/T iso ) : 



(3) 



where p' z is the boosted longitudinal momentum. 
One can again introduce auxiliary fields W 

Sf a (x;p) = -gW«(T,x\ m <t>,y)d? p) h{p^p v ) (4) 

that obey v ■ D W a (r 1 x l 7 r]; <j), y)\$, y = v^F a p, but now 

v a =p a /\p±\ = (cosh(y — rj), cos0, sin</>, T _1 sinh(?/ — 77)). (5) 

The induced current in the non-Abelian Maxwell equations is given by 

2-7T OO 

±D a (TF^)=.f(T,x\r,) = - HM[p™} f£j dyvPWiT^mhy) (6) 

-00 

with W = «¥j - Jr- sinh(y -rj)W n . 

In numerical simulations one now ends up discretizing velocity/rapidity space that has 
the form of a cylinder parametrized by </>, y instead of a sphere v 2 = 1. These numerical 
simulations have to start at a finite proper time r = To, both for practical reasons and 
for the physical reason that a plasma description makes sense at the earliest around r ~ 
Qj , where Q s is the saturation scale of perturbative QCD. Having chosen a background 
distribution function of the form ([3]), the only free parameters are Ti so , the (possibly 
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Fig. 1. Results 1191 from a 1D+3V real-time lattice simulation of non-Abelian plasma instabilities in 
Bjorkcn expansion, seeded by small initial rapidity fluctuations with a spectrum modelled after Ref. [13) . 
We show the proper-time dependence of the total chromo-field energy density and its individual compo- 
nents £ = £b t + £e t + £b l + £e l = £t + £l as well as the chromo-field energy gain rate R defined 
by R = dS/dr + 2£ t /t. 

fictitious) time when /o is momentarily isotropic, and the Debye mass following from /o 
at that point in time. 

In the results displayed in Fig. 1, t; so = O.lr, i.e. the distribution starts out already 
in oblate form and becomes increasingly oblate as time goes on. The mass parame- 
ter is matched to results from the saturation (colour-glass) scenario as explained in 
the appendix. In this numerical calculation, only the (most unstable) modes which are 
transversely constant and thus effectively 1+1-dimensional have been considered. These 
instabilities have been seeded by small rapidity fluctuations with a spectrum modelled 
after Ref. [13j . Given the findings in stationary anisotropic plasmas, the growth of non- 
Abelian plasma instabilities found here probably gives an upper bound on more generic 
cases. Considering that for Q s ~ 1 GeV for RHIC and ~ 3 GeV for the LHC and that 
thus the maximal time in Fig. 1 corresponds to some 20 fm/c for RHIC and 7 fm/c for 
the LHC, one finds an uncomfortably long delay for the onset of plasma instabilities at 
least for RHIC. Further studies of more generic initial conditions (including then also 
strong initial fields) is work in progress. 



Appendix A. Parameters of the numerical simulation 



For fixing the dimensionful parameters of the numerical simulation in a way that makes 
contact with heavy-ion physics, the so-called Colour-Glass-Condensate framework |14ll5j 
is invoked, and as starting time for the plasma phase r ~ Q~ l is chosen, where Q s is 
the so-called saturation scale. 

The only other dimensionful parameter, the Debye mass mu at the (fictitious because 
pre-plasma) time r; so , is fixed by assuming a squashed Bose-Einstein distribution function 
for the hard particle distribution function ([3]) with /i BO (p) = M(2N g )/{e p / T — 1) where 
N g = N% — 1 is the number of gluons and Af a normalization that is adjusted such that 
at t = tq the hard-gluon density of CGC estimates is matched. Since the expansion is by 
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assumption purely longitudinal, T is a constant transverse temperature, and it has indeed 
been found in CGC calculations that the gluon distribution is approximately thermal in 
the transverse directions, with T = Q s /d and dr 1 ~ 0.47 according to Ref. [15]. The 
normalization Af can then be fixed by following Ref. [16J . who write the initial hard- 
gluon density as 

N O 3 
Air z N c a s (Q s T Q ) 

where c is the gluon liberation factor, for which different estimates can be extracted from 
the literature. 

We adopted the value c = 2 In 2 w 1.386 from an analytical estimate in Ref. [UJ which 
turned out to be fairly close to the most recent numerical result c ~ 1.1 in Ref. |18) . With 
Tj so remaining a free parameter which determines how anisotropic the gluon distribution 
is at to, the normalization M is now fixed by n(ro)^ 2 - = n(r; so ) = 2< ^P AfN g T 3 . For a 
purely gluonic plasma, the isotropic Debye mass is given by m 2 D {Ti so ) = M^koc s N c T /3, 
which together leads to 



^(Tiso)r 2 (g s r )- 1 = ^ — « 1.285^-, N « 0.4628 a^Q.To)- 1 — 
6C(o) r iso r iso 7j 



ISO 



when c = 2 In 2 and AT C = 3. We adopt this value for our simulations where N c = 2, since 
in previous studies of the stationary anisotropic situation little difference was found 
between the SU(2) and the SU(3) case provided mo was the same |S]. With our choice 
of an initial anisotropy given by ro/ri SO = 10, equating To = Qj 1 and using units where 
r = 1, the above result corresponds to the value mn = 3.585 employed in Fig. Q] This 
corresponds to an energy density of hard particles £(t 1so = 0.1r ) = 7V87r 2 Qgd~ 4 /15 and 
(using Eq. (15) of [31]) £(t ) » 0.0789 £(t 1so ) w 0.0938 Q^QsTo)^^^ 1 . 
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